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Abstract. We present a method to search for large angular-scale correlations, termed topological signatures, in the angular 
distribution of cosmic objects, which does not depend on cosmological models or parameters and is based only on the angular 
coordinates of the objects. In order to explore Cosmic Microwave Background temperature fluctuations data, we applied this 
method to simulated distributions of objects in thin spherical shells located in three different multiply-connected Euclidean 3- 
spaces (T 3 , T n , and Ge), and found that the topological signatures due to these topologies can be revealed even if their intensities 
are small. We show how to detect such signatures for the cases of full-sky and partial-sky distributions of objects. This method 
can also be applied to other ensembles of cosmic objects, like galaxies or quasars, in order to reveal possible angular-scale 
correlations in their distributions. 
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1. Introduction 

Homogeneity property refers to the spatial distribution of cos- 
mic objects in the Universe, and to test it one needs to know 
the angular positions of these objects in the celestial sphere and 
their distances to us as well. Unfortunately, cosmological dis- 
tances are not directly measured, as they are calculated through 
the redshift-distance relationship, which involves cosmological 
parameters still not known with enough accuracy. So, precise 
distances to extragalactic objects are not currently available. 
Isotropy, instead, deals with the angular distribution of objects 
in the sky. Then, only two angular coordinates for each ob- 
ject are needed to determine whether or not they are isotropi- 
cally distributed around us. Nowadays such coordinates are ac- 
curately measured and there are available catalogs for several 
types of objects, like quasars, gamma-ray bursts, and the (sky 
pixels corresponding to the) Cosmic Microwave Background 
Radiation (CMBR) temperature fluctuations. In analyzing such 
catalogs, some questions arise: How can we recognize whether 
a given class of cosmic objects follows an exact or a quasi- 
exact isotropic distribution? How can we quantify small devia- 
tions from a purely isotropic distribution? How can we distin- 
guish a possible cosmological pattern or topological signature 
imprinted in these catalogs from pure statistical fluctuations? 
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The study of the large-scale (i.e. global) homogeneity and 
isotropy of the Universe encompasses its topology as well, 
since topological properties are global properties of 3-spaces. 
Multiply-connectedness is a topological property that tesselates 
a simply-connected space (like the Euclidean space K 3 ) pro- 
ducing multiple images of a given object. This property gener- 
ates distance and angular correlations (i.e. anisotropies) in the 
distribution of objects. The Small Universe hypothesis (Ellis & 
Schreiber 1986 1, which assumes that the Universe is a compact 
(i.e. multiply-connected with finite volume) 3-space that lies 
inside -at least partially- the Last Scattering Surface (LSS), 
has not been discarded by recent CMBR data. In view of this 
situation, we are motivated to study the anisotropic signatures 
of topological origin in the angular distribution of cosmic ob- 
jects, and for this we consider simulated ensembles of objects 
in different Euclidean multiply-connected 3-spaces. 

In order to search for these signatures in astronomical data 
sets, we developed a method based on histograms of angu- 
lar separations between pairs of cosmic objects. We call it 
the Pair Angular Separation Histogram (PASH) method. We 
are particularly interested in studying angular correlations in 
simulated catalogs of cosmic objects that share the same at- 
tributes as the photons of the CMBR temperature fluctuations, 
i.e., point-like objects located in the (comoving 3 -dimensional) 
thin shell representing the decoupling era, but viewed from 
Earth as located on the celestial sphere. As we shall show, 
the correlations originated by the isometry properties of the 
compact 3-spaces can be revealed by the method we propose 
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here. This method resembles the Cosmic Crystallography (CC) 
method (Lachieze-Rey & Luminet ll995l Lehoucq et al. 1996 
Fagundes & Gausmann I1998al I1998bl Gomero et al. I2000al 
20Q0b[ |Mni |M)2>. However, there is a basic difference be- 
tween them in that here we just look for angular correlations 
instead of searching for distance correlations as the CC method 
does. Our method is therefore independent of any cosmologi- 
cal model or parameters, while the CC method depends on this 
information to calculate the radial distance of cosmic objects 
from their measured redshifts. 

In Sect. 2 we present the basic properties of the 3-spaces 
representing our Universe in the Friedmann-Lemaitre cosmo- 
logical models (Friedmann 1922, Lemaitre ll927> . In Sect. 3, 
we present the geometrical and statistical method we developed 
to reveal correlations between pairs of objects distributed in a 
2-dimensional spherical surface. The treatment of this problem 
initiates with the theoretical description of a perfectly isotropic 
distribution of objects in a spherical surface S 2 . For this type of 
data we construct the normalized expected probability density 
for two objects be separated by a given angle, separation that of 
course ranges in the interval [Q,n]. This probability gives rise 
to the Expected Pair Angular Separation Histogram (EPASH) 1 . 
We also present the EPASH for partial-sky catalogs containing 
objects distributed only in polar-caps. In Sect. 4, we present the 
theoretical deduction of the topological signatures due to pure 
translational isometries, also known as Clifford translations, 
appearing in the PASHs. We illustrate these results, in Sect. 5, 
with numerical simulations using catalogs that contain multiple 
images of cosmic objects generated by the isometries of the fol- 
lowing Euclidean compact (i.e. multiply-connected with finite 
volume) orientable 3-spaces: T 3 , T„, and G(, (denoted by Q\, 
Qi, and Q(„ respectively, in Wolf's classification Wolf 1984 1. 
In this section we also show that performing the mean of sev- 
eral PASHs (or Mean Pair Angular Separation Histogram - 
MPASH) reduces significantly the statistical noise, which al- 
lows us to reveal the anisotropies due to a non-translational 
isometry whose tiny signature is present in the angular distri- 
bution of cosmic objects. We analized this problem considering 
MPASHs for full-sky catalogs as well as for polar-cap catalogs. 
Finally, in Sect. 6, we discuss our results and also the problem 
of finding suitable catalogs of cosmic objects from the available 
astronomical data. 



2. The geometry of the 3-spaces 

The manifold T 3 , best known as the three-torus, has its 
Fundamental Polyhedron (FP) usually -but not necesarily- rep- 
resented by an equal-sided cube, where its three pairs of faces 
are identified by pure (or Clifford) translation isometries. The 
manifold T„ has two pairs of faces identified by pure trans- 
lation isometries and one pair of faces identified after a rota- 
tion of 180°, the so-called screw-motion isometry. The man- 
ifold G(, has no pair of faces of its FP identified by a pure 



1 This result is independent of the three-dimensional geometry 
since all three isotropic geometries (i.e. spherical S 3 , Euclidean ft 3 , 
and hyperbolic H 3 ) possess spherical surfaces <S 2 as hypersurfaces. 



translation 2 . In analyzing the distributions of objects in these 
3-spaces we shall be considering the effect of not only pure 
translational isometries but also translation-with-rotation isom- 
etry, or screw-motion. For completeness, we find interesting to 
briefly review here the basic properties of compact 3-spaces of 
constant curvature. 

The large-scale structure of the Universe is described by the 
Friedmann-Lemaitre cosmological models. In these models, 
the local geometry satisfies the local homogeneity and isotropy 
properties and is described by the Robertson- Walker (RW) line 
element ds 2 = dt 2 — a 2 (t)dcr 2 , where a(t) is the scale factor, and 
do 2 = dx 2 + F k 2 (x) [d6 2 + sin 2 6d(p 2 ] is the line element of the 
t - constant spatial sections. The functions Fk(x) = sin^, x, 
sinh^ , where k is the normalized curvature parameter, which 
can take the values +1,0,-1, corresponding to the spherical, 
Euclidean, and hyperbolic geometry, respectively. It is known 
that the three-dimensional spatial sections (hereafter denoted 
by Mk) are manifolds of constant curvature obtained from the 
quotient Mk = Mk/T (see e.g. Wolf 1984 1, where T is a dis- 
crete group of isometries of the manifold Mk, without fixed 
points. The manifold Mk is termed the universal covering of 
Mk- Therefore, while Mk is a simply-connected manifold, Mk 
is a multiply-connected one. We would like to stress that the 
values k = +1, 0, -1 in the RW line element determines only 
the local geometry of Mk (equivalently that of Mk)- 

Since the simply-connected Mk and the multiply- 
connected Mk manifolds share the assumed local homogeneity 
and isotropy properties, we conclude that both manifolds are 
good solutions of the RW local geometry, i.e., the Einstein- 
Hilbert equations. Consequently, General Relativity (or any 
other metrical theory) can not distinguish between Mk and Mk 
because they look like each other locally. 

In a Universe with a multiply-connected 3-space Mk, and 
whenever the horizon scale is greater than half of the smallest 
closed geodesic of Mk, multiple images of a given cosmic ob- 
ject might exist due to the isometries of the group T. 3 These 
images originate distance correlations between the mapped 
cosmic objects and also produce angular correlations between 
those distance-correlated objects. As we shall see, these angu- 
lar correlations can be revealed by the PASH method described 
in the next section. 

3. PASHs 

3.1. EPASHs for full-sky catalogs 

In essence, a PASH is a normalized plot of the number of pairs 
of objects versus the angular distance between them. A formal 
construction of a PASH is as follows: let Sr c Mk be a ball of 
radius R centered at the origin of coordinates O, and involving 
spherical surface Sr, and containing an ensemble of (cosmic) 
objects. We call Sr the observable Universe, O represents the 
Earth, R_ is the radius of the LSS, Sr is the celestial sphere 

2 A picture of this FP can be seen, for instance, in Bernui et al. 
1 1998 1, where it is denoted by T4. 

Except when T = Identity, in which case Mk = Mk is a simply- 
connected manifold. 
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with radius R = R + — R + AR, where AR is the thickness of 
the decoupling era. Since we are interested in analyzing CMB 
data, our simulations consider ensembles of objects located in a 
thin shell of thickness AR, which are after projected onto the ce- 
lestial sphere Sr ignoring their radial coordinate. However, we 
point out that this method can be applied to any spatial distri- 
bution of cosmic objects. Given a set of well-defined selection 
rules 4 , a catalog C is a list of all objects in Sr which satisfies a 
sub-set of the whole set of rules. 

We start the description of the method by firstly showing 
how to obtain the PASH. Let's assume that we have an ensem- 
ble of objects that are, by assumption, isotropically distributed 
in Sr (see Eq.[5}. A full-sky map catalog C is a list of N ob- 
jects spreaded out in Sr with their corresponding angular coor- 
dinates. We divide the interval (0, n] in m bins of equal length 
6y - n/m, where each sub-interval has the form = (y. - , 
y. + y ] , i — 1, 2, . . . , m , with center in y. = (i — |) 6y . Now, 
we denote by <f>(y) the number of pairs of objects in C separated 
by a distance y e (0, n]. Thus, 



<*>(%■) 



N(N- \)5y 



(1) 



is the normalized counting of the number of pairs of ob- 
jects, separated by an angle y., that lies in the sub-interval 
Actually, Y!f=\ &(y-)8y = 1 , so we can average several func- 
tions <E>(y.), each one built from catalogs having a comparable 
number of objects. 

It may occur that some angular correlations of small inten- 
sity are present in a given PASH, but are not clearly seen. To 
reveal them, and since PASH and EPASH are normalized, we 
subtract the PASH from the EPASH in order to enhance possi- 
ble angular correlations. 

Let us now obtain the EPASH for the simply-connected 
manifold Mk- Let p,q € Sr c Mk be an arbitrary pair of 
objects listed in C with coordinates (r p , 8 P , ip p ), (r q , 6 q , <p q ), re- 
spectively. Then, the probability density that two completely 
uncorrelated objects 5 p and q be separated by an angle y e 
(0,7r],is 



V(y) = d*r p d*r q p(r p )p(r q ) S( A(p, q)-y), 



(2) 



where r p , r q are the vectors from O to the position of the ob- 
jects p and q, respectively. In Mo = 7? 3 , the angular separation 
between p and q is given by 

A(p, q) = arccos[cos 6 p cos 6 q + sin 9 p sin 6 q cos(<p p - <p q )], (3) 

where 6 p , 6 q e [0, 7r] and <p p , tp q e [0,2^]. Note that p(F p )d 3 r p 
is the probability of the object p to be in the volume cPr p . The 
probability of finding p in the whole Sr should of course be 
equal to 1, i.e., 



p(r p ) d 3 r p = 1 



(4) 



4 i.e. physical properties that are common to the cosmic objects in 
a given catalog, such as luminosity threshold, redshift, etc. 

5 this situation corresponds to the monopole term (/ = 0) in a mul- 
tipole decomposition. 
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Fig. 1. EPASH (smooth curve) and PASH (noisy curve) for a 
full-sky simulated catalog with ,/V = 200 objects in Sr; m = 
400 bins. 

Thus, a purely isotropic and normalized density distribution of 
objects observed in Sr reads 



P(r p ) 



S(r p -R) 



pJ 4nR2 
Finally, integrating we obtain 

= j sin y . 
Hence, the EPASH is simply defined by 

<&expec,ed(y,) = T~ f Pfy) • 



(5) 



(6) 



(7) 



However, if the interval 7, is small enough (as shall be consid- 
ered here), a suitable approximation for the EPASH is 



^expected(7i) ~ ^(y,) • 



(8) 



We also note that P(y) defined in Eq. satisfies the normal- 
ization property, J Q f P{y)dy = 1, which is a useful condition 
since it let us to perform the mean of an arbitrary number of his- 
tograms obtained with catalogs containing a comparable num- 
ber of objects. 

In Figs, n an d ID we show the EPASH together with the 
corresponding PASHs obtained from numerical simulations for 
full-sky catalogs with ,/V = 200 and ,/V = 2 x 10 3 objects, re- 
spectively. 

3.2. EPASHs for polar-cap catalogs 

Our galaxy poses a serious problem to the use of full-sky 
CMBR maps and for other sky catalogs as well, because it con- 
taminates part of the celestial sphere around the Galactic plane. 
This fact makes somewhat useless the approach we developed 
studying the angular distribution of cosmic objects listed in 
full-sky catalogs. For this reason, we studied also the PASHs 
for catalogs containing objects located in polar-caps (i.e. above 
or below the equator of a spherical distribution). 
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Fig. 2. EPASH (smooth curve) and PASH (noisy curve) for a 
full-sky simulated catalog with N = 2 x 10 3 objects in Sr\ m 
400 bins. 

Consider a catalog with objects located in a spherical cap 
surface S pc c S 2 , centered, for instance, around the z-axis, 
which is assumed to be aligned with the North Pole. We call 
this ensemble a polar-cap catalog Cpc- In this case, any object 
with coordinates 



x — R cos if sin 
y = R sin^j sin6>, 
z — R cos 0, 



(9) 



where <p e [0, 2n] and the azimuthal angle 8 e [0, Gq], with 8q < 
n/2, is in the surface Spc- Thus, the isotropic and normalized 
density distribution of objects in S pc is 

&(8 o -0 p )6(r p -R) 



P(r p ) = 



(10) 



IPC 



where Ape = 2nR 2 {\ - cos#o) is the area of the polar-cap. 
Using this information in the definition of an EPASH (Eq. |3 
and after the integration of the 5 Dirac-function, we obtain 

Ppciy, Q ) =-^- f sin B p de p f iin q d0 q I(0 p , q ; y), 
A pc Jo Jo 



(ID 



where 



pin r*2n 

I(0 P , Q q \ r) = I I d<p p d<p q d(A(p, q) - 
Jo Jo 



y). 



It can be shown that this integral is equal to 

1(8 8 ) 47rsiny 

p ' 9 ' sin 8 P sin ft, si\\a(8 p ,8 q \y) ' 

where the function a is defined by 

cos y - cos 6 p cos 8 q 



o{8 p , 8 q ; y) = arccos 



sin Op sin 8 q 



(12) 



(13) 



(14) 



and lies in the interval [0, n]. Due to this condition, P and q 
should satisfy the relation 



cos y - cos P cos q 
sin 8 P sin 8 q 



e [-1,1]. 



(15) 



Fig. 3. EPASHpc (smooth curve) and PASH (noisy curve) for a 
polar-cap simulated catalog with 0q = 70°, Npc = 332 objects 
(10 3 objects in the full-sky simulated catalog, m = 400 bins); 
cr = 0.02799. 

Equivalently, this constraint for the values of P and q can be 
written as the following two inequalities: 



p +0 q > y, 



(16) 



\0 p -0 q \ < min{0 o ,r}. 



Thus, the probability density to observe two objects in a polar 
cap (defined by the angle 0q) with an angular separation y is 

4;rfl 4 sin yf e ° f e ° d0 q 
Ppc(% Oo) = ~ 2 d0 p . q , (17) 

A 2 pc Jo 'J sma(O p ,e q ;y) 

with the conditions given in Eq. (I16> . The double integral in 
Eq. O can be numerically solved, but it is convenient to do 
it by parts considering carefully the conditions imposed by the 
inequalities in Eq. d!6l >. Thus, given 8q, we can obtain the prob- 
ability density Vpc as a function of the angular separation y. 

Recently, the analytical expression for the expected proba- 
bility density that an arbritary pair of objects in S pc be sepa- 
rated by y degrees, i.e. EPASH for polar-caps (EPASHpc), has 



been found by Teixeira ( 2003 ), who obtained 
4-nR 4 sin y 



A 2 

A PC 



■ ( arcsin[ {(cos 8q + cos y) 



-\/cosy - cos 20o } / {sin#o (1 + cos 8q) y/l + cosy } ] 
+ (1 -2 cos ft ) arccos [cot 8q tan(y/2)]). 



(18) 



To illustrate this result, we simulate two full-sky catalogs with 
= 10 3 and N - 4- x 10 3 objects, respectively. We show 
in Figs.|3]and|4]the EPASHpc together with the corresponding 
PASHs obtained from these simulations, where the number of 
objects in the polar-caps were Npc = 332 and Npc = 1337, 
respectively. 

3.3. MPASHs for full-sky and polar-cap catalogs 

Suppose that K comparable catalogs Ck, k — 1,2, K, of ob- 
jects located in the celestial sphere are available for analysis. 
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Fig. 4. EPASHpc (smooth curve) and PASH (noisy curve) for 
a polar-cap simulated catalog with 80 = 70°, Npc = 1337 
objects (with 4 x 10 3 objects in the full-sky simulated catalog, 
m = 400 bins); cr = 0.00661. 



Fig. 5. Difference between PASH and EPASH plotted in Fig.Q] 
for a full-sky simulated catalog (K — 1) with N = 200 objects 
and m = 400 bins; cr = 0.04418. 



For a given m (common to all the PASHs), we calculate the 
functions 



**(7i) = 



N k (N k -l)6y 



(19) 



yejj 



where N k is the number of objects in C k , and <f> k (y) is the num- 
ber of pairs of objects in C k with angular separation y. Then, 
the MPASH is defined by 



1 K 

<<">(y s )> = °M 
K ti 



(20) 



Indeed, we know that performing the average of Pair 
Separation Histograms in the CC method, the result is that the 
statistical fluctuations are reduced by a factor proportional to 
1 / but the topological information is preserved (Gomero et 
al. I2000al feOOObl 12001 1. In the next section we shall see that 
for the PASHs a similar situation occurs, that is, the greater 
the number of PASHs considered in the average, the better the 
approximation of the MPASH to the corresponding EPASH. 

3.4. Noise in PASHs 

In Figs.|5]and[6]we show the difference between the PASH and 
the EPASH, already plotted in Figs. 1 and 2, respectively. The 
noise of a PASH, as observed in Figs.|5]and[6] is simply given 
by 



rj(y.) = 0(y.) - <$> expected {y?) , 



(21) 



and is obtained from the very definition of the EPASH. In fact, 
given a catalog C with N objects, we can quantify the statistical 
noise in a PASH by: 



1 

\ m ^ 



(22) 



P- 0.005 



1.0 1.5 2.0 2.5 

Angular Separation (radians) 



Fig. 6. Difference between PASH and EPASH plotted in Fig.EI 
for a full-sky simulated catalog (K - 1) with ,/V = 2 x 10 3 
objects and m — 400 bins; cr = 0.00453. Note that the scale 
here is different from that of Fig.|5] 



For the cases shown in Figs.|5]and[6] we evaluate the statistical 
noise using Eq. (|22ll obtaning cr = 0.04418 for N = 200, and 
cr = 0.00453 for N = 2 x 10 3 respectively. 

In Fig. we show the difference between the MPASH 
and the EPASH for a numerical simulation of K — 10 4 full- 
sky simulated catalogs in 7? 3 , with a mean number of (AO = 
1 /KYjf = i Ni = 200 objects per catalog. The statistical noise is 
cr = 0.00045. Fig.[8]shows the difference between a MPASH, 
obtained from a numerical simulation of K - 10 4 full-sky 
catalogs with (N) = 20 objects each, and the correspond- 
ing EPASH; in this case cr = 0.00448. In Fig. |5] we show 
the difference between a MPASH and the EPASH, where the 
MPASH was obtained from a numerical simulation of K = 150 
full-sky simulated catalogs with (N) = 100 objects each, but 
with m = 600. In this case cr = 0.00926. 
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Fig. 7. Difference between MPASH (for a numerical simula- 
tion of K — 10 4 full-sky simulated catalogs in H 3 , with (N) = 
200 objects and m = 400 bins) and EPASH; cr = 0.00045. 
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-0.04 I i i i , i , i i , i , i i i i , , i i i , , , , i , , , , i ,_ 

0.0 0,5 1.0 1.5 2.0 2.5 3.0 

Angular Separation (radians) 

Fig. 9. Difference between MPASH (for a numerical simula- 
tion of K — 150 full-sky simulated catalogs in fl 3 , with (AO = 
100 objects and m = 600 bins) and EPASH; cr = 0.00926. 
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Fig. 8. Difference between MPASH (for a numerical simula- 
tion of K — 10 4 full-sky simulated catalogs in ft 3 , with (N) = 
20 objects and m = 400 bins) and EPASH; cr = 0.00448. 

Thus, when we compare the cases shown in Figs. [6] and [8] 
(equal m and cr, different values for N and K) with the cases 
shown in Figs.|5]and0(equal m and N, different values for cr 
and K), we finally arrive at 

0-45 ^ .... 
cr - — . (23) 

(AO 

We can verify this relationship with the aid of the case shown 
in Fig.|5] in which m = 600. In fact, using Eq. d23l . we obtain 
cr = 0.009, while calculating the right hand side of Eq. ( 1221 
with the values of <E>(y,), i = 1, ...,600, obtained from Fig. [5] 
results cr = 0.00926. 

In Figs. |3 and [6] as well as in Figs. [3] and 0] we notice the 
absence of correlations between pairs of objects for the exactly 
isotropic distribution in ft 3 generated using random numbers 
for the coordinates. Thus, it is clear from these simulations that 
as K — > oo and/or N — » oo the statistical noise measured by 



cr — > 0. The MPASH is related to the statistical noise and this 
fact is quantified in Eq. fl23i . Although the statistical noise is 
inherent to any numerical simulation, it can be controlled in 
three ways: i) by obtaining catalogs with a large number of 
elements, in order to produce PASHs with a small noise, since 
the statistical noise depends on the number of objects; ii) by 
computing the MPASH to significantly reduce the statistical 
noise; iii) by choosing a suitable number of bins m. 

4. Topological signatures in PASHs 

In any multiply-connected 3-space Mo whose FP lies inside the 
ball Sr, there might exist images of a given object in a shell 
of thickness AR = R - R , whenever the horizon scale R is 
greater than half of the length of the smallest closed geodesic 
of Mo (Lachieze-Rey & Luminet [T995l Lehoucq et al. IT996I 
Fagundes & Gausmann [T998al I1998bt Gomero et al. I2000al 
2000b 2001 2002 1. These images give extra contributions to 
the PASH due to the isometries of Mo- 

By definition, the probability density for two objects p, q 
in C be correlated by a translational isometry gi (we represent 
this correlation by q — gi(p) ) is 

p^i)= V JL f f d\ p d\ q p(r p )p{r q )x 

6{d{p,q)-i)5(r q - gi (r p )), (24) 

where I is the distance variable for objects in Br, I e (0,27?]; 
d(p,q) = \f p - r q \ is the 3-dimensional Euclidean distance be- 
tween p and q, and the coefficients v,- = v„. are the number of 
images generated by the isometry g jm Moreover, n is the density 
number of objects inside the volume under investigation (in the 
present case a thin shell), and Af is the number of objects in this 
volume. 

Consider the pure translational isometries in Mo = T 3 , 
with a FP represented by a cube of side L. We know that the 
distances between pairs correlated by Clifford translations are 
independent of the location of the object (see Gomero et al. 
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20021, i.e., V/7 e S R , the distance \gip)\ = d(p,gi(p)) = A t 
is independent of p. Thus, for this manifold, the translational 
isometries g,, i = 1,2,... are such that for A\ — L we obtain 
V! = 6, for A2 = V2i we obtain V2 - 12, for /1 3 = s/3L we 
obtain v 3 = 8, etc. 

Let's now apply the above definition to the topological im- 
ages correlated by the Clifford translation isometries g, and lo- 
cated in a thin shell of thickness AR. In this case, 



P(f p ) 



&(R-r p )®(r p -R-) 
f[R3-Rl] 



(25) 



Using this expression in Eq. ( 1241 . and after some calculations, 
we arrive at 



(26) 



which, due to the 6 Dirac-function, appears as an acute peak 
localized at I — A t in any pair distance-separation histogram, 
and where 



Pi = 



3 v t AR (R + R-Y 
8 Aj N (R 2 + RR + Rl) 



(27) 



gives the height of the peak. We notice that the probability for 
these isometries to show up in pair separation histograms is 
proportional to the thickness AR of the shell and inversely pro- 
portional to N, the number of objects in the shell. 

In Euclidean multiply-connected 3-spaces with Clifford 
translations, the objects located in a thin shell and correlated in 
distance by the Clifford translational isometries gj (of the cov- 
ering group F) give rise to a definite signature in the PASHs. 
In fact, for shells of small thickness (like that corresponding to 
the decoupling era where AR/R =s 0.003), these angular cor- 
relations between pairs appear in the MPASHs as small peaks 
and are approximately located at the angular scales 



2 arcsin 



(28) 



whenever R > A mi „/2, which is the requirement for the exis- 
tence of topological images in the observable universe. Thus, 
it is not difficult to verify that the minimum and the maximum 
angle {y 1 "'" and yf ax , respectively) subtended by two objects 
located in a thin shell and separated by a isometry distance At 
are, respectively: 



sin(yf /2)-A. 



sin(y" M 72): 



2R- 



(29) 
(30) 



We observe that the probability for angular correlation between 
objects is a direct consequence of the existence of the N(N - 
1 )/2 pairs located in the thin shell and correlated by the Clifford 
translation Aj with the probability given in Eq. ( 1271 . In fact, the 
probability density that a pair of objects be correlated by to the 
isometry gi at an angle y, is 



'(r) = Pi 



0(r _ r ™» )0(r »««_ r ) 



(31) 



where Ay ? = y"™ x - y" m - Thus, the topological signatures 
coming from the isometries gi contribute to the EPASH (i.e. 
Eq.[7} as 



O ft . ,, see 

expected 7 



expected 



SyJ 



by 



(32) 



Therefore, the angular correlations due to Clifford translational 
isometries appear in the PASHs as triangular peaks with basis 
Ay and height Pi/Sy. 

The effects of the non-translational isometries can also be 
calculated using Eq. J24i . However, in such a case, the inte- 
grals are cumbersome to handle, so we prefer to show their ef- 
fect in MPASHs with plentiful numerical simulations. For this 
reason, we decided to work with the G(, manifold, which does 
not present Clifford translations in the histograms whenever 
R/L < V2/2, L being the arista of the cube representing the 
FP. 

5. Numerical simulations 

We performed numerical simulations aiming to reveal topolog- 
ical signatures due to the isometries of the compact flat man- 
ifods r 3 , T„ and Ge for full sky and for polar-cap simulated 
catalogs. 

It is difficult sometimes to recognize a topological signa- 
ture in a single PASH due to the presence of statistical noise. 
As a matter of fact, we found that topological signatures due to 
Clifford translations appear in PASHs as small triangular peaks 
whose heights are inversely proportional to the number of ob- 
jects in a given catalog (see Eq. 1261 . We saw that it is con- 
venient to divide a catalog of N objects in K sub-catalogs of 
(N) — N /K objects, where the objects in a sub-catalog share 
common physical properties, and then perform the MPASH of 
the K resulting PASHs. By doing this, we magnify in two ways 
the presence of the topological signatures: reducing the num- 
ber of objects per catalog, N — > (N), the height of the peak 
increases by a factor N/ (N) = K; and obtaining K comparable 
catalogs we can perform the MPASH to reduce the noise by a 
factor 1 / V^- The net result is to increase the signal-to-noise 
ratio (SNR) of the topological signatures by a factor V^. This 
explains why the strategy of dividing in sub-catalogs the origi- 
nal catalog with a large N, and performing the MPASH works 
better than considering just the PASH of the original catalog. 
For 3D catalogs, that is, when the 3 spatial coordinates of the 
objects are provided, the Collecting Correlated Pairs method 
(Uzan et al. 1999) is suitable to enhance the SNR for the signa- 
tures originated by the holonomies of the 3-space Mk- 

In the simulations presented in this section we assume that 
the FP is a cube of side L — 1, and the origin of coordi- 
nates is in the center of this cube. Moreover, for T 3 and T n , 
we considers = R + = 1, R = 0.997, i.e. AR = 0.003. The 
case of Gf, is doubly special: firstly, in order to avoid the pres- 
ence of Clifford translations appearing at A — V2/2, we as- 
sume R = 0.7; secondly, because the non-translational isome- 
tries have a very small signature in the MPASH, we assume 
R = 0.69, i.e. AR = 0.01. 

In order to normalize and compare the results appearing in 
the different MPASHs shown below, we fix the values of m, 
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and consequently that of 6y, in all of them. Thus, we consider 
that the interval [0, n] is divided in m = 400 bins, each one of 
width 6y = n/m = 0.00785. 

In general, to clearly reveal the signatures in the MPASH, a 
good SNR is of the order of 10; in particular, regarding Clifford 
translational isometries, this implies that the number of cata- 
logs should be such that 

^ ~ 10 o\ (33) 
by 

Thus, using the expression for P,, Eq. i27\ . and the numer- 
ical values for the parameters involved, we obtain that, for 
K i 10 4 , the topological signatures due to Clifford translations 
are clearly revealed in the MPASHs. 

5.1. MPASHs for full-sky maps 

In this section, we show how to generate full-sky catalogs in 
order to construct the MPASHs to reveal topological signatures 
in the angular distribution of cosmic objects. We assume in all 
the simulations presented here that the FP is a cube of side L = 
1, centered at the origin of coordinates. To construct a PASH 
through a numerical simulation, we first randomly generate a 
set of three Euclidean coordinates x,y,z e [-0.5, 0.5] for N see d s 
objects located inside the FP. 

In order to be sure that the distribution generated in this 
way is uniform, we test the randomness of the simulation using 
the Pair Separation Histograms for the simply-connected case 
7? 3 , according to the theoretical relationship expected for such a 
case given by Bernui & Teixeira ( 1999 1. After that, these N see d s 
objects are mapped to the whole ball Sr via the translational 
isometries gj —>/},= 1, V2, V3, ... of T 3 , giving rise to Nj 
objects (including the N see ds objects) in Sr. Then, we select the 
N objects located inside the thin shell, i.e. those objects with 
radial distance satisfying R < r\ < R, i = 1, ...,N. We project 
them on the involving sphere Sr and calculate all the angular 
distances between the N(N - l)/2 pairs of objects. 

We then count how many of these angular distances are 
in the interval i = 1,2,..., 400 . Finally, we normalize the 
PASH and plot the number of distances versus the angular dis- 
tance y. For the simulations showed in Figs. [TO] and where 
we consider the case Mq = T 3 , we use R — 1, R = 0.997. 
In this case, due to the ratio between the volumes of Sr 
(= 47r/3) and T 3 (= 1) we should obtain images correspond- 
ing to the first three Clifford translations: A\ = 1, Aj = V2, 
and A3 = V3. As a consequence of this, the MPASH presents 
angular correlations that are approximately located at y™" ^ 
1.047, 1.571, 2.094 radians, respectively, with the correspond- 
ing peak heights Pj/dy ~ 1.14xl0~ 2 , 1.6xl0~ 2 , and0.9xl0~ 2 . 

For the simulation showed in Fig.^] where we study the 
case Mo = T„, we also used R = 1, R = 0.997. However, in 
this case, the Clifford translations contributing to the topologi- 
cal angular correlations studied in the previous section are just 
two: A\ = 1, A2 = V2, which generate only the images v\ = 
4, V2 = 4, respectively. According to this, the heights of the 2 
peaks appearing in Fig.[Q]are PJby ~ 7.6 x 10~ 3 , 5.4 x 10" 3 , 
and are approximately located at y™'" ^ 1.047, 1.571 radians, 
for i = 1,2 respectively. In fact, as we can observe in Fig. 1131 



1 1 




Angular Separation (radians) 

Fig. 10. MPASH together with EPASH for K = 10 4 full-sky 
simulated catalogs in T 3 (the FP is a cube of side 1); (AO - 100 
objects. 
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Fig. 11. Difference between MPASH and EPASH correspond- 
ing to Fig.fTUI 

there is no peak at y^'" ^ 2.094 because there is no transla- 
tional isometry at A — V3 in T n . 

In Fig. Elwe present the MPASH and the EPASH for the 
case of the 3-space Mq = G(,, and in Fig. ^] we show the cor- 
responding difference between them. In this case, we assumed 
again that the FP is a cube of side 1, R = 0.7, R = 0.69, i.e. 
AR = 0.01, (N) - 20 objects, and K = 4 x 10 5 . 

Comparing Figs. II II and 1131 with Fig. [H] we observe that 
the heights of the peaks are much smaller in the case of screw- 
motion isometries than in the case of Clifford translation isome- 
tries. However, our simulations show that with a suitable num- 
ber of catalogs, the screw-motions, as well as the pure transla- 
tional isometries, are revealed in the MPASHs as small peaks. 

5.2. MPASHs for polar-cap maps 

We present numerical simulations in which the objects are lo- 
cated in a thin spherical shell of thickness AR; the radial dis- 
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1.0 1.5 2.0 

Angular Separation (radians) 



Fig. 12. MPASH together with EPASH for K = 2 x 10 4 full-sky 
simulated catalogs in T„ (the FP is a cube of side 1), (N) - 100 
objects. 




Fig. 13. Difference between MPASH and EPASH correspond- 
ing to Fig.fT2l 




Fig. 14. MPASH together with EPASH for K = 4 x 10 5 full-sky 
simulated catalogs with objects in (the FP is a cube of side 
1), R = 0.7, R = 0.69 (i.e. AR = 0.01), (N) =* 20 objects. 
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Fig. 15. Difference between MPASH and EPASH correspond- 
ing to Fig. 1141 



tance of the object to the origin of coordinates is R < r < R , 
&o = 70°. We performed these calculations for both T 3 and T„ 
3-spaces. 

A noticeable fact, shown in Figs.[^[El andl 1 5 Icorresnond- 
ing to the numerical simulations done in the precedent sub- 
section, is that full-sky catalogs produce small distinguishable 
peaks in the MPASHs when the isometries are Clifford trans- 
lations (which appears only in Euclidean and spherical geome- 
tries; for details see, e.g., Gomero et al. 2000a] I2000bl |2"0"0"T1 
2002 1, and even smaller signatures for the non-translational 
isometries. 

For completeness we perform now numerical simulations 
in order to test whether these results are valid for partial-sky 
catalogs, namely for polar-cap catalogs, and we do this for 
the cases when the 3-space is T 3 and T n . Observing Figs. El 
and^] we notice that the topological information does not dis- 
appear when we consider a suitable polar-cap catalog instead of 
a full-sky catalog. 

Regarding our simulations, we remark that the values as- 
sumed for R and L means that the observable Universe includes 
a small number of multiple images, generated by the isometries 
of the manifold, for each seed-object present in the FP. This as- 
sumption is within the limits imposed by COBE to the scale of 
the Small Universe considering a T 3 manifold (Stevens et al. 
[19931 de Oliveira-Costa & Smoot fT995l . 



6. Discussion and conclusions 

We presented a method to study the angular distribution of cos- 
mic objects. It has the great advantage of being independent of 
cosmological models or parameters. Through the analyses of 
the PASHs, we investigated the effect that different topologi- 
cal properties, termed Euclidean isometries, produce in the an- 
gular distributions of objects. We performed these analyses in 
full-sky as well as in polar-cap catalogs. 

We concentrate on the analysis of objects located in a thin 
spherical shell, considering only their angular positions in the 
celestial sphere in order to explore CMBR data. We have con- 
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Fig. 16. MPASH together with EPASH for a polar-cap catalog 
in T\ with O = 70°, (N) - 40 objects, K = 4 x 10 4 . 



,'Vi 1 
|i 'i i 



i .0 1 .5 2.0 2.5 

Angular Separation (radians) 



Fig. 17. Difference between MPASH and EPASH plotted in 
Fig. [16] 



sidered four different cases regarding the 3-space where the 
shell is embedded: the simply-connected "R 3 , and the three 
topologically different multiply-connected fiat 3-spaces: T 3 , 
T„, and Ge- 

In the simply-connected case, topological angular corre- 
lations are clearly absent. This is evinced by plotting the 
difference between the MPASH and the EPASH, where only 
statistical fluctuations appear, as seen in Figs.[5]to[9] indepen- 
dently of the number of simulated catalogs used to produce the 
MPASH. For the multiply-connected 3-spaces analyzed, an- 
gular correlations due to the isometries involved in each spe- 
cific case indeed appear. Thus, the topological imprints cor- 
responding to translational isometries manifest themselves as 
small and independent peaks, i.e., one peak for each isometry. 
On the other hand, the signatures corresponding to the non- 
translational isometries (which is the case for the Gf, manifold) 
appear as a less intense but extended angular correlation: from 
y™" =i 1.571 radians to y™ m 2.094, as can be seen in Fig. [131 



Fig. 18. Difference between MPASH and EPASH obtained for 
a polar-cap catalog (with #o = 70°) for the case where the 3- 
space is T„; (N) - 40 objects and K = 4 x 10 4 . The axis of 
symmetry of the polar-cap coincides with the axis of the screw- 
motion isometry. 

It is also important to notice that the topological informa- 
tion of multiply connected 3-spaces does not disappear when 
we consider a suitable polar-cap catalog instead of a full-sky 
catalog, that is, provided that such polar-cap is sufficiently 
large and conveniently oriented in the sky to map multiple im- 
ages. This fact is clearly observed in Figs. I17landll8l This is 
a very useful result because in general astronomical data are 
presented in catalogs covering only partial regions of the sky. 
In order to improve the SNR and reveal small topological im- 
prints in the angular distribution of cosmic objects, the strategy 
is to divide the original catalog in several sub-catalogs and then 
perform the MPASH. It could be also useful to divide the full- 
sky data in a set of antipodal polar-caps and compare them to 
similar sky patches in a statistically isotropic Universe (Bernui 
et al. l2004l 

We have shown that a topological signature, no matter how 
small it is, can always be revealed whenever one achieves a 
suitable SNR in the difference between the MPASH and the 
EPASH. As shown in our simulations, the small topological 
signatures of T 3 and T„ are clearly revealed when averaging 
K = 4 x 10 4 catalogs with (N) - 40 objects each, from an orig- 
inal catalog of 1.6 x 10 6 cosmic objects. For G^, the number 
of cosmic objects should be five times larger. For the case of 
CMBR maps which are always contaminated by Galactic fore- 
grounds, the polar-cap analysis described in this paper would 
be useful to reveal possible angular correlations in these data 
sets. 

Acknowledgements. We thank A.F.F. Teixeira, G.I. Gomero, M.J. 
Reboucas, and C.A. Wuensche for many fruitful discussions. A.B. 
thanks CNPq (PCI/DTI7B fellowship) and TWAS for the partial finan- 
cial support, and the hospitality of the Centro Brasileiro de Pesquisas 
Ffsicas visited under the TWAS-Associateship Scheme at Centres of 
Excellence in the South. TV. was partially supported by CNPq grant 
302266/88-7-FA and FAPESP grant 00/06770-2. We also thank the 
anonymous referee for the valuable comments on the manuscript. 



Bernui & Villela: Topological signatures in the angular distribution of cosmic objects 



11 



References 

Bernui, A., Gomero, G.I., Reboucas, M.J., & Teixeira, A.F.F. 1998, 

Phys. Rev. D, 57, 4699 

Bernui, A., & Teixeira, A.F.F. 1999, astro-ph/9904180 
Bernui, A., Villela, T, & Ferreira, I. 2004, Int. J. Mod. Phys. D, 13, 
1189 

de Oliveira-Costa, A., & Smoot, G.F 1995, ApJ, 448, 477 
Ellis, G.F.R., & Schreiber, G. 1986, Phys. Lett. A, 115, 97 
Fagundes, H.V., & Gausmann, E. 1998, |astro-ph/9811368| 
Fagundes, H.V., & Gausmann, E. 1998, Phys. Lett. A, 238, 235 
Friedmann, A. 1922, Z. Phys., 10, 377 

Gomero, G.I., Reboucas, M.J., & Teixeira, A.F.F. 2000, Phys. Lett. A, 
275, 355 

Gomero, G.I., Reboucas, M.J., & Teixeira, A.F.F. 2000, Int. J. Mod. 
Phys. D, 9, 687 

Gomero, G.I., Reboucas, M.J., & Teixeira, A.F.F. 2001, CQG, 18, 
1885 

Gomero, G.I., Teixeira, A.F.F, Reboucas, M.J., & Bernui, A. 2002, 

Int. J. Mod. Phys. D, 11, 869 
Lachieze-Rey, M., & Luminet, J.-P. 1995, Phys. Rep., 254, 135 
Lehoucq, R., Lachieze-Rey, M., & Luminet, J.-P. 1996, A&A, 313, 

339 

Lemaitre, G. 1927, Ann. Soc. Sci. Bruxelles, 47, 29 

Stevens, D., Scott, D., & Silk, J. 1993, Phys. Rev. Lett., 71, 20 

Teixeira, A.F.F. 2003, physics/0312013 

Uzan, J.-Ph., Lehoucq, R., & Luminet, J.-P. 1999, A&A, 351, 766 
Wolf, J. A. 1984, Spaces of Constant Curvature (Delaware: Publish or 
Perish Inc.) 



